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Abstract In this review we discuss some observational aspects and theoretical models 
«^ [ of astrophysical collisionless shocks in partly ionized plasma with the presence of non- 

thermal components. A specific feature of fast strong collisionless shocks is their ability 
to accelerate energetic particles that can modify the shock upstream flow and form the 
M^ ■ shock precursors. We discuss the effects of energetic particle acceleration and associated 

magnetic field amplification and decay in the extended shock precursors on the line 
and continuum multi-wavelength emission spectra of the shocks. Both Balmer-type 
^ ■ and radiative astrophysical shocks are discussed in connection to supernova remnants 

I , interacting with partially neutral clouds. Quantitative models described in the review 

Q ■ predict a number of observable line- like emission features that can be used to reveal the 

physical state of the matter in the shock precursors and the character of nonthermal 
^ . processes in the shocks. Implications of recent progress of gamma-ray observations of 

supernova remnants in molecular clouds are highlighted. 
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1 Introduction 



Plasma flows with coUisionless shocks are found in a number of energetic space ob- 
j ects, such as supernova remnants (SNRs) interacting w i th atomic o r molecular clouds 
( McKee and Hoilenbachlll980l:lDraine and McKe3ll993l: iHendboiol'). He rbig-Haro ob- 
jects ( Hartigan et al.Ml987l : iGiannini et al.M2008l: iTesileanu et al.M2009r ), winds from 
protostars and young s tellar objects ( MassLet al. 200g), recurrent novae ( Evans et alj 
120071 : iBode et al.ll2007l : iTatische ff and Hernanz 2007 ) , and also in accretion flows in the 
vicinities of galactic nuclei (IFarage et alJ 201G| ). 

Already in the early studies of iRavmondl ( 19791 ) and lShuU and McKed (|l979f l the 
ionization state of the pre-shock flow was recognized as an important feedback param- 
eter that influences the dynamics of the flow, and hence, t he spectrum of c ontinuum 
and line emission produced in the post-shock zone. While Raymond 1 19791 ) took the 
pre-shock ionization state as a free parameter, IShuU and McKed l| 19791 ) included the 
ionizing flux from the post-shock to obtain a pre-shock ionizati on state self-consistently 
with t he temperature and ionization proflles in the post-shock. iDopita and Sutherland! 
( 19961 ) furthered the calculation of the precursor photoionization, and presented emis- 
sion line spectra assuming an equilibrium with the ionizing flux from the shock. 

Recent observations of forward shocks in galactic SNRs indicate a su bstantial role 
of no nthermal components in the energy budgets of post-shock flows (e.g.,[Helder_et_alJ 
I2OI2I ). While the fraction of the energy dissipated in the shock itself may be smaller 
for slower radiative shocks, the enhanced radiative cooling in their post-shocks leads 
to stronger compression so that cosmic ray and magnetic pressure can dominate over 
thermal pressure in the zones where most of the observed optical and IR emission 
is produced. There is also a growing consensus on that in both parallel and per- 
pendicular shocks the diffusing cosmic rays may generate unstable plasma currents 
and turbulence in the pre-shock zone, thus allowing for efficient particle accelera- 
tio n, shock modiflcation, and substantial enhancement of the total compression ra- 



tio IIBlandford and Eichlen 



1987; 



__. Ijones a nd Ellisonl ll99ll : iMalkov and Drurvlliooil : iBelll 

12004 : iBvkov et al. II 20 121 ). I Raymond et al.. (,1988bf ) analyzed a shock in the Cygnus Loop 
SNR and found that the nonthermal pressure exceeds the thermal pressure by an or- 
der of magnitude in the zones where the [S II] lines are formed, though they were 
unable to distinguish between magnetic and cosmic ray con tributions. Neutral parti- 



cles could affect the p r ocesses of particle acceleration (e.g.. iDraine and McKee 199^ 



O'C Drurv et al.lll996': 'Malkov et al]|2005l: iBlasi et all [2013: 
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2012l : lHelder et al.ll20 12: Mor lino et al.ll2012al). magnetic fleld ampliflcation and plas ma 
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heating in the upstream region (e.g.. lBvkov and Toptvginll2005l : iReville et al.ll2007l ). 

In this paper we review physical mechanisms governing evolution of supersonic 
and superalfvenic flows with coUisionless shocks in partially ionized plasmas with non- 
thermal components. Section [5] describes the photoionization precursors that help de- 
termine the neutral fraction of the pre-shock gas. In SectionOHa signatures of interac- 
tion of neutral particles with fast shocks are depicted. In Section |4] we discuss specific 
plasma heating processes due to dissipation of CR-driven magnetic fluctuations in the 
shock precursor. A consistent pre-shock ionization structure is given in Section [S] with 
an aim to emphasize processes that critically influence optical and IR spectra emitted 
from the flows. The escape of accelerated energetic particles may affect the gas com- 
pression in the shock upstream and plasma temperature in the downstream resulting 
in formation of radiative shocks of velocities that may well exceed that of the adiabatic 
shocks. We discuss the line spectra of radiative shocks modifled by CR acceleration 



effect in Section [S] Specific features of CR acceleration in supernova remnants inter- 
acting with molecular clouds are discussed in Section |6] and their gamma-ray emission 
spectra are outlined in Section [T] A brief summary is given in Section [S] 



2 Photoionization Precursors 

The neutral fraction in the pre-shock gas often depends on photoionization by radiation 
from the shock itself. There are two regimes, non-radiative shocks are those in which 
gas is heated at the shock but does not have time to radiatively cool, while radiative 
shocks radiatively cool from the post-shock temperature, efficiently converting their 
thermal energy into radiation. Naturally, the latter produce more ionizing photons if 
the post-shock temperature exceeds about 150,000 K. The scale of the photoionization 
precursor is given by the phot oabsorpt ion cross section of hydrogen and the neutral 
density in the pre-shock gas, so it is typically 10 or 10 cm. In supernova remnants, 
shocks faster than 300 km/s are typically nonradiative, while slower shocks are generally 
radiative. 

Non-radiative shocks produce some ionizing photons because some of the atoms 
passing through the shock are excited before they are ionized in the hot post-shock 
plasma. In particular, at high temperatures each He atom produces a few photons in the 
He I 58.4 nm and He II 30.4 nm lines. Since He makes up about 10% of the gas, the pre- 
shock region will be by 30% to 40% ionized by these photons. EUV emission from the 
hot post-shock gas can further ionize the plasma. The photons are relatively energetic, 
so they heat the gas to around 16,0 00 K. Emission from such prec ursors outside Tych o's 
SNR and N132D were reported bv lChavamian et al.1 (|2000l ) and lMorse et all (|l996h . 



As mentioned above, and as described more fully in Section [5] radiative shocks 
faster than about 100 km/s produce large fluxes of ionizing photons, and those above 
about 120 km/s produce enough photons to fully ionize H in the pre-shock gas. 
However, it should be kept in mind that the equilibrium between pre-shock ioniza - 
tion state and ionizing flux f rom the shock, as assumed in IShull and McKed ( 19791 ): 



iDopita and Sutherland! ( 19961 ). does not always hold. In particular, an SNR shock can 



evolve more rapidly than the recombination time in the pre-shock gas, or the transverse 
scale of the shock may b e smaller than the ups tream photon mean free path, as it is 
in a Herbig-Haro object ( Ravmond et al.lll988al ). 



3 Interaction of Neutrals with Fast Shocks: Ha Signatures 

Neutral H atoms are not affected by the electromagnetic flelds or plasma turbulence in 
a coUisionless shock. They pass freely through the shock transition, whose thickness is 
given by the ion inertial length or the proton gyroradius, which are typically ~ 10 cm 
in the interstellar medium. However, they become ionized in a hot downstream plasma 
but, because of the low ISM densities, they travel ~ 10 -10 cm downstream be- 
fore it happens. This is smaller than electron-proton thermal equilibration scales, and 
therefore, the emission from neutrals may contain information about the post-shock 
conditions of electrons and ions before relaxation. 

The hydrogen atoms are ionized either by collisions with protons and electrons or by 
a charge exchange with protons. The latter process produces a population of neutrals 
with a velocity distribution similar to that of the post-shock protons and, when they are 
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Fig. 1 Ha profile of a 350 kin/s shock at the outer edge of the Cygnus Loop, obtained with 
the HECTOCHELLE spectrograph on the MMT. The narow component is 30 km/s wide 
(FWHM), and the broad component is 225 km/s wide. Full line is the observation, dashed line 
is the Ha from the filament, and dotted line is the sky background. The feature near 6569A is 
an airglow line. 



excited, they produce broad Lyman a nd Balmer emission lines whose velocity widths 
are comparable to the shock speed ( Chevalier and Raymond! 1 19781 : iRavmondl Il99ll : 
iHend uOlOl : [France et al.l uOllfl . Neutrals that are excited before the charge transfer 
produce a narrow component whose width is given by the pre-shock temperature (see 
Figure[T]). The broad-to-narrow intensity ratio indicates the electron temperature, while 
the broad component width directly measures the proton temperature. Thus one can 
determine the electron-ion temperature ratio at the shock (Ghavamian and Schwarz this 
volume). Consequently, the Ha profile provides both the post-shock tem perature and 
Te/Ty ratio. These quantities can be used to determine the shock speed pmith et al.l 
I1994I : iHester et aPl 19941 : IChavamian et al.ll200ll : Ivan Adelsberg et al.ll200a l. Under the 
assumption that only a small fraction of the shock energy goes into cosmic rays, the 
shock speed can be determined from the post-shock proton and electron temperature 
using momentum and energy conservation laws. However, if cosmic ray acceleration is 
efficient, then Te and Tp c ombined with a m odel of efficient cosmic ray acceleration 



can yield a shock velocity ( Helder et al.ll2009l ) though there is some ambiguity due to 
a contribution to the narrow component from a shock precursor (see Section [3. 2p . 



The Balmer line profiles can also indicate more exotic plasma processes. First, if the 
post-shock proton velocity distribution is not Maxwellian, it will leave an imprint on the 
broad component line profile. Second, the strength and width of the narrow component 
can indicate heating and excitation in a shock precursor. We discuss signatures of these 
departures from a simple shock picture. 



3.1 Non-Maxwellian Velocity Distributions 

II cosmic rays are accelerated in SNR shocks, the particle velocity distribution is mani- 
festly non-Maxwellian, taking the form of a Maxwellian core with a power law tail. Such 
distributions are conveniently parameterized with the k function (|Pierrard and Lazaij 
I2OI0I ). The K distribution approaches a Maxwellian as k approaches infinity, while the 
power law tail becomes harder and includes a larger fraction of the particles as k 
approaches its lower limit of 1.5. 

A different sort of departure from a Maxwellian can arise because of the neutral 
atoms themselves. When a neutral becomes ionized, it is moving at 3/4 the shock 
speed with respect to the downstream plasma and magnetic field (assuming a shock 
compression factor of 4). It therefore behaves like a pickup ion in the solar wind 
( Gloeckler and Geisd I2OOII ) . Initially, all these newly-formed protons have the same 



velocity with respect to plasma in which they are immersed. This velocity separates 
into a gyration speed around the magnetic field plus a component along the field, and 
since all the particles have the same initial velocity (assuming that the shock speed 
is large compared to the pre-shock thermal speed), they form a ring beam in velocity 
space. That distribution is unstable, and it rapidly evolves into a bispherical shell in 
velocity space, a lens- like shape that depends o n the ratio of the Alfven speed to the 
initial particle speed ( Williams and Zanklll994h . The protons in the bispherical shell 



can then experience charge transfer to produce observable neutrals. The process of 
ionization, interaction with the magnetic field, and subsequent re-neutralization is the 
same sequence of events that produce the "IBEX ribbon" seen in H atoms from beyond 

the heliopause (|McComas et akluOOgf ). 

Th e Ha profiles produced by the pickup ion process were computed bv lRavmond et al.l 
( 20081 ) . The neutral fraction in the pre-shock gas must be substantial in order to pro- 



duce the Balmer line filaments seen at the edges of some supernova remnants and 
pulsar-wind neb ulae, and therefore the neutrals can have a significant effect on the 
shock structure. iHeng et al.l (|2007l ') computed the gradual transiti on over the charge 
transfer length scale, but found little effect on the overall dvnamics. lOhira et al.l ( 2003 ) 



examined the streaming instabilities created by the relative motion of the post-shock 
plasma and the protons created downstream (this effect is stro nger in parallel shocks ) 
and predicted substantial amplification of the magnetic field. IRavmond et al.l (|2008l ) 
calculated the wave energy produced as the ring beam relaxes to a bispherical distribu- 
tion, and showed that it could significantly affect the elect ron temperature if the waves 
couple efficiently to electrons. lOhira and Takaharal ( 20ld ) considered the modification 



of the shock structure, in particular weakening of the subshock, and the effects on the 
cosmic ray spectrum. 

The observational evidence for any departure of the proton distribution from a 
Maxwellian is still ambiguous. For nearly all the observed Ha broad components Gaus- 
sian fits are adequate, but in general the lines are so faint that even the determination of 
the line width is subject to large statistical uncertainty. A very deep exposure of a bright 
knot in Tycho's SNR yielded the only profile that shows a clear non-Gaussian broad 
component, but there are several possible interpretations ( Ravmond et al.ll20ld ). The 



departure could be attributed to a power law tail, however, the slope of this tail is very 
hard, which indicates extremely efficient particle acceleration. The profile could also be 
matched by the combination of a pickup ion component with an ordinary Maxwellian 
component, though this requires a fairly high pre-shock neutral fraction. A third pos- 
sibility is that some of the neutrals in the shock precursor aquire a kinetic temperature 



about 1/2 of the post-shock temperature. Clearly, this implies a precursor to be both 
hot and thick. Finally, one could produce the observed profile by adding contributions 
from separate shocks along the line of sight. However, this possibility would require 
quite a large velocity difference, and therefore a significant density contrast between 
the two regions. At higher shock spee ds, the velocity dependen c e of the charge tr ansfer 
cross section can distort the profile ( Heng and McCravl 120071 : iHeng et al.l 120071 1. but 



that is unlikely to account for the Tycho observations. 



3.2 Shock Precursors 

The precursor of a collisionless shock wave is a region upstream of the shock transition 
in which the plasma conditions (velocity, density, temperature, magnetic fields, ioniza- 
tion state) are affected by photons or superthermal particles streaming ahead of the 
shock front. 

Shock wave precursors can be produced by ionizing photons emerging from the 
downstream region or by broad component neutrals that leak back through the shock. 
Additionally, in shocks that produce accelerated super-thermal charged particles, pre- 
cursors can be produced as a consequence of this acceleration process. We first consider 
the observational evidence for precursors, then the physics of the three mechanisms and 
the implications for system parameters. 

The first evidence for precursors to SNR shocks came from the observation that 
the widths of the Ha narrow components in several shocks were unexpectedly large. 
ISmith et al.l ([199J) measured nar row component wid ths of 25-58 km/s in four LMC 
Balmer-dominated remnants, and iHester et al.l ( 1994 ) found a 28-35 km/s width in a 



Cygnus Loop shock. These line widths correspond to temperatures of about (2-7) x 
10 K, and i n stati c equi Ubrium, hydr o gen is fully ionized at those temperatures. Thus 
Smith et al.l (|l994l l and iHester et al.l (|l994l ') concluded that the gas must be heated 
in a narrow region ahead of the shock. The region must be thick enough so that 
charge transfer can heat the neutrals, but thin enough that the neutrals do not become 
i onized. Subsequent obse rvations have revealed line widths of 44 km/s in Tycho's SNR 
( Ghavamian et al.l I2OO0I ') and 30-42 km/s in RCW 86 and Kepler's SNR, with only 



SN 1006 showing a narrow compone nt width compatible w it h a temperature of 10,000 



K and a significant neutral fraction ( SoUerman et al.ll2003f l. lNikolic et al.l ( 20131 ) have 



presented Integral Field Unit spectra of a section of the Hq filament in SN 1006, 
and they have found that a precursor makes a substantial contribution to the narrow 
component of Ha. 

The precursor thicknesses must be about l" in order for the neutrals to be heated , 
and they have b een spatially resolved in the Cygnus Loop ( Fesen and Itohl 1 19851 ; 
IHester et al1ll994l ') and in Tycho's SNR JLee et al.ll2007l . l2010l ). 



3.3 Physical Interpretation of Shock Precursors 

Three ideas have been put forward to explain observed shock precursors, a photoioniza- 
tion precursor, a cosmic ray precursor, or a precursor created by fast neutrals leaking 
from the post-shock region out ahead of the shock. 

The photoionization precursor heats the electrons, and it has a very large length 
scale. In non-radiative shocks, each He atom passing through the shock produces about 



2 He II A304 photons and 1 to 2 He I A584 photons. These photons can ionize about 
30% of the pre-shock hydrogen, and because the photons are relatively energetic they 
deposit considerable heat in th e electrons. Such a prec ursor has been reported ahead 
of Tycho's supernova remnant ( Ghavamian et al.ll2000l ). However, because they heat 



the electrons and extend over large length scales, photoionization precursors do not 
explain the observations discussed in Section [3.21 

A precursor is an integral part of diffusive shock acceleration models, in which 
an MHD t urbulence scatters energetic particles moving away from the shock back 



towards it ( Blandford and Eichleilll987l l. Its thickness is given by kcr/Vs, where Kcfl 
is the cosmic ray diffusion coefRcient. While k.(jr is of order 10 cm /s in the ISM, it 
must be somewhere closer to the Bohm limit near the shock in order for particles to 
reach high energies, perhaps 10 cm /s. That implies a precursor thickness of order 
10 cm, or about l" for nearby SNRs. The heating in such a precursor can occur by 
dissipation of the turbulence itself. It will be especially strong if neutrals are present 
to provide ion-neutral wave damping ^ which can limit the wave intensity and therefor e 
the highest particle energies reached ( Draine and McKeelll993l : IO'C Drurv et al.lll996l l. 



The heating will also be strongest in shocks that accelerate cosmic rays efficiently. If the 
precursor heating is strong, it can reduce the M ach number of the subshock, changing 
the cosmic ray spectru m and pressure (see, e.g.. lVladimirov et al.ll2008l : IWagner et al.l 
l2009l : IVinket al.ll2O10l ). 

Several theoretical models have considered t he effects of neutra l -ion c ollisions on 
the dynamics and temperature of the precursor. iBoulares and Coxl ( 19881 ) computed 



the heating and i onization in cosmic r a y pre cursors for the modest shock speeds in 
the Cygnus Loop. lOhira and Takaharal (|2010r i considered neutrals interacting in the 
precursor, treating them as pickup ions. The pickup ions can strongly heat the plasma 
and affect the jump conditions at the subshock, and because of thei r high velocities 
they a re preferentially injected into the diffusive acceleration process. iRavmond et al.l 
(|2011h computed the temperature and density structures of precursors and the resulting 
Ha emission. If cosmic ray acceleration is efficient, the Ha profiles are strongly mod- 
ified, both in broadening the narrow component and decreasing the broad-to-narrow 
intensity ratio. The Ha profiles reported so far are consistent with moderate ratios of 
cosmic ray pressure to gas pressure, but not with ratios above about 0.4. Most recently, 
iMorlino et al.l ( 2012b|) have included a parameterized form of wave damping along with 



wave growth due to the cosmic ray pressure gradient for a more self-consistent treat- 
ment of the cosmic ray spectrum in the presence of neutrals and the resulting Ha 
profiles. They also find that an interme diate velocity width com ponent can arise, sim- 
ilar to the complex profiles predicted in lRavmond et al.l ( 201H ). 



The other suggested mechanism for producing a precursor is leakage of some of 
the broad component neutrals out through the shock front. These neutrals have the 
thermal speeds of post-shock protons and a bulk speed of V5/4 away from the shock. 
A significant f raction of them, of order 10%, can ove rtake the shock, and of course they 
pass through ( Smith et al.lll994l : lHester et al.lll994l l. They carry a substantial amount 



of energy, but it is unclear how much of that energy heats the precursor plasma. 
iLim and Ragal 1 1996h found that they became ionized to form a beam of fast protons. 



but that the protons are swept back t hrough the shock without transferring much en- 
ergy to the bulk plasma. lOhiral (|2012l 'l constructed four-fluid models and showed that 
the le aking partic l es dec elerate the upstream flow and affect the subshock compression 
ratio. iBlasi et al.l ( 20121 ) computed the energy and momentum exchange b etween neu- 



trals and ions on both sides of the subshock, using a procedure like that of lHeng et al.l 



(120071) to compute the neutral velocity distribution and the flux of neutrals from the 
post-shock region to the precursor. Then they used the resulting precursor structure 
to compute spectral slopes of accelerated particles. The authors accounted for neutrals 
ionization by protons only, though phot o ionization processes can affect the ionization 
states for some shock param eter space. iMorlino et al.l (2012a|) used the kinetic simu- 



lations of iBlasi et al.l (|2012h to calculate the Ha emission produced in the precursor. 
Its intensity depends on the poorly constrained efficiency of electron-ion equilibration, 
since the heat is deposited in the ions, but the electrons excite the line. They find that 
the profiles can deviate strongly from Gaussian and that in some circumstances an 
intermediate width comp onent can arise . They suggest that this could explain profiles 



observed in Tycho's SNR. IMorlino et al.l (2012b|) have extended this work by including 



the cosmic rays along with the neutrals that escape upstream through the shock. 



4 CR precursor heating and the post-shock temperature 

Models of coUisionless shocks with large sonic and Alfvenic Mach numbers {A4s » 1 
and Alfvenic A4a. » 1) show that, through the first-order Fermi acceleration mech- 
anism, a small minority of particles could gain a disproportionate share of the energy 
and populate the high energy tail of particle distribution. The energetic particles can 
penetrate far into the shock upstream gas, to create an extended shock precursor illus- 
trated in Figure [2] The cold gas in the shock precursor is decelerated and pre- heated 
by fiuctuating magnetic field dissipation on a scale that is about c/ugh times larger 
than a mean free path of an energetic particle A*. Shocks in coUisionless supersonic 
flows produce a complex multi-scale structure of the relaxation region with an extended 
precursor and sub-shock of a modest sonic Mach number Algub ~ 3. 

This section discusses theoretical models of strong coUisionless shocks and their 
implications for the observations of shock precursors discussed in Section 13.21 



4.1 The structure of a non- linear shock precursor in the presence of self-generated 
MHD turbulence 

Consider a strong, plane-parallel coUisionless shock in a plasma of finite /3 = TMaZ-'^s 
(/? is the ratio of plasma pressure to magnetic pressure). The distribution function 
of nonthermal particles and the bulk flow profile in the shock upstream region are 
sensitive to both the total upstream compression ratio rtot and the subshock Mach 
number Msuh- Direct numerical simulations of the CR- modified shock by particle- 
in-cell technique are nonfeasible by now because of the wide dynamical range of the 
simulation that requires extreme computing resources. Nevertheless , an approximate 
iterative approach (e.g., within the Monte Carlo mo del discussed in 
(200S) or semi-analytical kinetic models developed bv lMalkovl (19971); 
( 



Vladimirov et al 



Amato and Blasi 



20061 ') 1 can be used to derive the steady-state distribution function consistent with 
the shock compression. These approximate models assume some diffusion model and 
parameterize the microphysical processes of magnetic field amplification and plasma 
heating. The exact calculation of the CR escape flux Qosc that determines the total 
upstream compression ratio rtot can be p erformed only in fully no nlinear self-consistent 



simulations. The Monte Carlo model of I Vladimirov et al.l ( 20081 ) describes the escape 



of particles from the shock with an assumed free escape boundary far upstream of the 
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Fig. 2 A schematic view of a CR-modified shock propating to the right. 



shock. The distance to the free escape boundary is a free parameter of tire simulation 
that controls the maxim um energy of accelerated particles and the escaping CR flux 
(see Ellison et al.l (| 19961 1 for more details of this method) . 

The role of accelerated particles in the shock precursor depends on the poorly 
known mechanisms with which energetic particles can transfer energy to the ther- 
mal pool. Howev e r, in order to numerically estimate the impact of such a process, 
IVladimirov et al.l (|2008l l parameterize the rate of plasma heating upstream. In order 
to do that, they assume that CR particles produce, through the resonant CR stream- 
ing instability, strong fluctuations of magnetic fields upstream of the shock. Another 
assumption in the model is that these fluctuations are immediately dissipated into the 
heat of the thermal gas, and that the rate of turbulence dissipation is proportional 
to the rate of turbulence generation. The dimensionless parameter Qff in the model 
i s the ratio of dissipation to generation rate. Figure [3] adopted from I Vladimirov et al.l 
(120081) shows the self-consistently calculated flow speed, effective magnetic field and the 
temperature in the precursor of a strong shock with turbulence dissipation for several 
values of a^ ■ This work demonstrated that even a small rate of turbulence dissipation 
can significantly increase the gas temperature in the precursor, and thus increase the 
rate of injection of thermal particles. At the same time, the spectrum of high energy 
particles is only modestly affected by dissipation in these simulations. 

In a follow-up work of lVladimirov et al.l ( 2009! ) the effect of the microphysics of the 
fluctuating magnetic fields in the shock precursor was modeled in a more physically 
realistic manner. That work used the balance equation of the turbulence spectral energy 
density 

^+uVrVy + Vkn = 5(a.,k)-Lt„b- (1) 

Here I/turb(2;i k) represents the dissipation of turbulence, the turbulence energy injec- 
tion rate is S{x, k) = ^{x, \i.)W{x, k, i), and 7(2;, k) is t he rate of wave energy amplifica- 
tion by the CR-driven instabilities (see for a discussion iBvkov et al.ll2012l : ISchure et al.l 
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Fi g. 3 The profile o f a CR -modificd shock simulated with Monte-Carlo nonlinear DSA model 
by I Vladimirov et al.1 ll2008l ). The solid, dashed and dotted lines correspond, respectively, to a 
fraction, ajj = 0, 0.5 and 1.0, of the magnetic energy that goes into heating of the thermal 
gas. The plotted quantities are the bulk flow speed u{x), the effective amplified magnetic field 
Bgf{(x) and the thermal gas temperature T{x). The shock is located at a; = 0. 



bold ). The fastest short- wavelength CR-driven instability studied by iBelll (|2004h has 
the following growth rate for the wavevector along the mean magnetic field: 



7CR = 2«Afc|| 



for l/^gi < k < kc, 



(2) 



where vj^{x) = Bq / y^ in p{x) is the Alfven speed, c is the speed of hght, Bq is the 
far upstream magnetic field directed towards the shock normal, p{x) is the thermal 
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plasma mass density, rgi{x) is the gyroradius of the least energetic current gener- 
ating CR, the critical wavenumber kc{x) = 47rj(i(a;)/(c_Bo), and the local diffusive 
electric current of CRs responsible for the instability, jd{x), is determined via a Monte 
Carlo simulation. The turbulence energy density flux n(r, k, t) describes anisotropic 
cascading, i.e., the transfer of turbulence energy from long to short wavelengths (Kol- 
mogorov type cascade) as well as the inverse cascade (see, e.a.. lMonin and Iaglomlll975l : 
IZhou and Matthaeus|[l990l : IVermalbOOJ : IZhoull20ld ). 

Modeling the balance, spectral distribution, and dissipation of turbulence energy 
with equation 11} allows us to construct more physically realistic models of shock 
precursors and account for the Balmer line observations discussed in Section [S] This 
equation must be supplemented by theoretical models of the growth, dissipation and 
spectral flux of MHD turbulence produced by CRs. Sections 14.21 and 14.31 discuss some 
of the ongoing work in this area. 



4.2 MHD turbulence damping in partially ionized plasmas 

Ion-neutral collisions may dominate the frictional damping of strong magnetic oscilla- 
tions i n a cold photo-ionized pla s ma (T ~ 2 x 10 K) of a precursor of a CR-modified 



shock ( Draine and McKed Il993l : lO'C Drurv et al. igQa). In a cold partia lly ionized 
plasma the generalized Ohm's law (see, e.g.. lBraginskiilll965l : ICowlindll976h 



E+-uxB = J- + — j X B + ^ ^'° B X (j X B). (3) 

c a ruec nimiC^ 

results in enhanced effective magnetic diffusion v^t in the induction equation 

47rcr_J 
where 

1 1 f'^b'^ 

— T = - + 7- (5) 

Here F is the mass fraction of the neutrals, t^^ is the mean time between the ion-neutral 
collisions. The neutrals in cold magnetized plasma of supernova shock precursors may 
strongly affect both the growth and the dissipa t ion of CR-driven waves (see, e.g., 
iBvkov and ToptvginI boOSi : iMarcowith et al]|2006l : iReville et al.llJOOTl '). In a precursor 
of a quasi-parallel shock the CR-driven turbulence is basically incompressible and the 
energy dissipation rate f™ in the low /3 limit can be estimated as 

F = aj^ fc|| , 

while the ion-neutral dissipation of the compressible fluctuations t hat are associated 
with long-wavelen gth CR-driven oblique wave instabilities (see, e.g.. lBvkov et al.ll201ll : 
ISchure et al.ll2012l ) are determined by 

F"' = alf k'^. 

The magnetic turbulence dissipation term in Eq. fl} due to effective magnetic 
diffusivity (Joule dissipation) is determined by 

-^tuib(2:, k) = u^fk'^W{x, k, t). (6) 
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To characterize a fraction of the magnetic energy dissipated in the precursor of a 
CR-modified shock it is instructive to estimate the minimal wavenumber km of the 
fluctuations to be strongly dissipated by ion-neutral friction while advecting through 
the CR-shock precursor of the scale length Lcif ~ c/i'sh ^* ^s it is illustrated in 
Figure [2] where 

t, k1y,h£R^l, (7) 

Then the maximal wavelength to be damped is 

Am = 2-K/km ~ 2tiFM^^ X ^CTjaA*. (8) 

For a strong supernova shock of Ma ~ 100 and F ~ 0.1 one may get Am ^10 cm. 
Note, however, that the length Am is comparable to the ion-neutral collision length that 
is ~ [F Tii ain\~ , if the ion number density in the precursor Ui <1 cm~ . Depending 
on the spectral shape of the CR-driven magnetic turbulence in the shock precursor, the 
dissipated fraction ranges from a few percent for a Kolmogorov-type spectrum to about 
10% for the flatter spectra of strong turbulence. We neglected here the MHD wave 
dissipati on due to the th ermal conduction and viscosity assuming a cold plasma case 
(see, e.e:.. lBrae:inskiilll965l . for a thorough discussion). The magnetic power dissipated by 
the ion-neutral collisions in the low /3 plasma mainly heats ions. The coUisional damping 
discussed above results in a true irreversible conversion of the magnetic turbulence free 
energy into the thermal plasma energy. To the contrary, the coUisionless turbulence 
dissipation processes that we are going to discuss below can heat electrons, thus being 
not completely irreversible. The coUisionless heating process just increases the energy 
of qu asi-thermal plasma components, but in general some coUisionality (e.g., iRamod 
[201l|) is required to increase the entropy and to reach the equilibrium distributions of 
plasma species. 



4.3 CoUisionless heating of ions and electrons by magnetic turbulence 

The CR-driven turbulence source S{x,\<i) in the shock pr ecursor is expected to be 
anisotropic in k. The fast CR-current driven instability of iBelll (200J) as well as the 



long- wavelength instability of iBvkov et al.l ( 201l|) in the Bohm diffusion regime are 



mainly amplifying the modes with t he wavevectors along the local magnetic field. The 
acoustic instability of IPrurv and Fallg (jl98^) driven by the CR-pressure gradient is 
also anisotropic. The strong anisotropy of the magnetic turbulence is observed in the 
solar wind where the outward fiux significantly exceeds the ingoing one. 

The solar wind is one of the best laboratories to stud y anisotropic magnetic turbu- 
lent dissipation and coUisionless plasma heating (see, e.g., lLeamon et al.|l998|; Sahraoui et al.l 
I2OO9I : IPetrosvan et aLlboid : lAlexandrova et al1l201lh . Recentlv. Isahraoui et al.l (|2009l ') 
reported Cluster spacecraft measurement providing direct evidence that the dissipation 
range of magnetofiuid turbulence in the solar wind extends down to the electron scales. 
Namely, they found two distinct breakpoints in the magnetic spectrum at frequencies 
fp — 0.4 Hz and fe = 35 Hz, which correspond, respectively, to the Doppler-shifted 
proton and electron gyroscales. Below fp, the spectrum follows a Kolmogorov-type 
scaling of power-law index about -1.62. For fp < f < fe a, second inertial range with 
a scaling index -2.3 was established. Above fe the spectrum has a steeper power law 
-4.1 down to the noise level of the Cluster detectors. The authors advocated a good 
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agreement of the results with theoretica l predictions of a quasi- two-dimensional cas- 
cade into Kinetic Alfven Waves fKAWV lChen et alJ (|2012l l presented a measurement 
of the scale-dependent, three-dimensional structure of the magnetic field fluctuations 
in the inertial range of the solar wind turbulence. The Alfven -type fluctuations are 
three-dimensionally anisotropic, with the sense of this anisotropy varying from large to 
small scales. At the outer scale, the magnetic field correlations are longest in the local 
fluctuation direction. At the proton gyroscale, they are longest along the local mean 
fleld direction and shortest in the direction perpendicular to the local mean field and 
the local field fiuctuation. The compressive fiuctuations are highly elongated along the 
local mean field direction, although axially symmetric in the perpendicular direction. 
Their large anisotropy may explain why they are not heavily damped in the solar wind 
by the Landau damping. 

Anisotropic wavenumber spectra that are broader in wavenumber perpendicular 
to the background magnetic field than in the parallel one are expected even in the 
symmetric MHD case studied by iGoldreich and Sridhail (| 19951 ) where the oppositely 
directed waves carry equal energy fluxes . The main features of t he th e symmetric (but 
anisotropic) MHD turbulence model bv IGoldreich and Sridhail ( 19951 ') are: 



(i) a critical balance between the linear wave mode periods and their nonlinear 
turnover timescales in the inertial-interval energy spectrum; 

(ii) the 'eddies' are elongated in the direction of the field on small spatial scales 

1/3 2/3 
with the scaling fcii oc fcp kJ , where fco is the wavenumber corresponding to the outer 

scale of the turbulence. 



The three-dimension al simulations by Cho and LazarianI ( 2004 ) revealed the basic 



features of the IGoldreic h and Sridhail ( 1995l l model in the electron magnetohydrody 



namic turbulence. The kinetic Alfven wave and whistler fluctuations are likely to make 
an important co n tribution to the turb u lence below the proton g y roscale (see, e.g., 
Garv et all I2OI2I : ISaito and Garvl I2OI2I : iBoldvrev and Perej I2OI2I : iMithaiwala et all 



2OI2I ). Particle-in-cell simulations show that the anisotropic whistler turbulence heats 



the electrons in the parallel direction as predicted by the linear theory and that in 
the low P plasmas the magnetic wavenumber spectrum becomes strongly anisotropic 
with spectral index in the perpendicular direction close to -4 (see, e.g., iGarv et all 
I2OI2I : ISaito and Garvl I2OI2I ) . Microscopic tw o-dimensional p a rticle- in-cell simulations 
of whistler turbulence were carried out by ISaito and Garvl (|2012l l in a collisionless 
homogeneous magnetized plasma to study the electron and ion heating dependence 
on the plasma magnetization parameter /3. They demonstrated that at higher values 
of P the magnetic energy cascade in the perpendicular direction becomes weaker and 
leads to more isotropic wavenumber spectra. The electron energy ratio between parallel 
and perpendicular components becomes closer to unity at higher /3. Three-dimensional 
particle-in-cell plasma simulations of decay of initial long wavelength whistler fiuctua- 
tions into a broadband, anisotropic , turbulent s pectr um at shorter wavelengths via a 
forward cascade were performed bv lGarv et all (|2012l ). The simulations demonstrated 
a picture qualitatively similar to that in 2D but somewhat stronger anisotropy of 
the resulting 3D turbulence comparing to the similar 2D runs. They showed a clear 
break in the perpendicular wavenumber spectra qualitatively similar to that mea- 
sured in the electron scale fluc tuations in the solar wind. Earlier lOuataertl (|l998l l and 
iQuataert and Gruzinovl (| 19991 ) discussed the beta-dependence of particle heating by 
turbulence in advection-dominated accretion flows. They found that for /3 ~ 1, i.e. ap- 
proximately equipartition magnetic fields, the turbulence primarily heats the electrons. 
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For weaker magnetic fields, the protons are primarily heated. The division between elec- 
tron and proton heating occurs between 5 < /3 < 100, depending on unknown details 
of how Alfven waves are converted into whistlers at the proton gyroscales. 

The cascade of Alfven waves, which are weakly damped down to the scale of the 
proton gyro-radius k \ p j ~ 1 is a subject of gyro-kinetic models (see for a review 
ISchekochihin et al.ll2009f l. The models can be of interest for the parallel shock precursor 
heating assuming an efficient cascading of the CR-driven magnetic fluctuations down 
to the fcj^pj; ~ 1 regime. The cascading still remains to be demonstrated since the CR- 
current driven modes are very different from the standard Alfven waves. If the cascading 
occurs, then the continuity equation ^ can be reduced to the equation for bf.{k±) = 
kj_ J dfcii W(k) - the energy density of the a nisotropic magnetic turbulence as a function 
of the perpendicular wavenumber (see, e.g.. lHowesl2010l : ICranmer and van Ballegooiienl 
I2OI2I I: 

^ + k^^ = sik^)-rik^)bl (9) 

where the energy injection rate is S (non-zero only at the driving scale kj_ — ko), the 
linear energy damping rate is F, and e(fc^) is the energy cascade rate. To specify the 
energy cascad e rate both advection and diffusion in the wavenumber space models are 
used (see, e.g.. lCranmer and van Ballegooiienl ( 20121 ) for a recent discussion). Assuming 



critical balance at all scales and using the energy cascade rate in the form of the 
advection in the wavenumber space: 

e{k^) = C;^/^kj_uibl 
iHowea (2010|) obtained the steady state solution for the energy cascade rate as 



where C'l and C2 are the dimensionless Kolmogo rov con s tants (Ci = 1.96 and C2 ~ 
1.09) and eq is the rate of energy input at /cq. iHowesI ( 201(3 ) used the normalized 



eigenfrequencies io{kj_) from the linear gyroki netic dispersion relat ions and the damping 
rates Fs due to different plasma species from lHowes et al.l ( 20061 ) and Eq. (|10p (where 
s = i, e), to calculate the spectrum of heating by species 

0,s(fcx) = cf 'C2(r,/zj)6(fc^)/fc^. 

The ion damping peaks at k^pi ~ 1, while the electron damping peaks at k^pi 2> 1 
unless Ti/Te < me/rrii. The energy that passes through the peak of the ion damping at 
k±Pi ~ 1 would lead to electron heating assuming both the cascading and the damping 
times are shorter than the advection time through the shock precursor. Then the total 
(integrated over kj_ pi > 1) ion-to-electron heating rate due to the kinetic dissipation 
of the turbulent cascade, Qi/Qe{l3i,Ti/Te) can be approximately fitted with 



+ /?f lm,F, _i/0. 



where ci = 0.92, ca = 1.6/{Fi/Fe), C3 = 18 -f Slog (T,/Te), and p = 2 - 0.21og(T,/Te). 
A slightl y better fit for Fi/Te < 1 occurs with the coefficients C2 = 1.2/(T;/re) and 
C3 — 18 ( HowesluOKj ). The model is valid for the parameter range 0.01 < Pi < 100 and 
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0.2 < Ti/Te < 100. The heating rate Qi/Qe is an approximately monotonic function 
of Pi that is only weakly dependent of Ti/Te- 

The simplified model of the electron and ion heating discussed above assumed 
the nonlinear coUisionless cascading from the energy containing scale of wavenumber 
k*Rg ~ 1 where the amplitude of amplified magnetic field is SB* that is determined 
by the gyroscale Rg of the energy containing accelerated particles to the thermal ion 
gyroscale kj_pi ~ 1. In the shock frame the advection time of the incoming plasma 
through the CR precursor r^dv ~ -^* c/i'sh- Therefore for efficient heating of the plasma 
species the cascading time Tc must be shorter than Tadv If the cascading time is 
determined by the turn-over time of the energy containing magnetic "vortex" t^ ~ 
fc*(5_B*/y47rp then the condition of efficient plasma heating in the shock precursor by 
the CR-driven turbulence can be written as 

'4^>1' (12) 

where rj = X*/Rg > 1. The amplitude of amplified magneti c field likely scales with the 
shock velocity as (5_B* ex u^j^, where 1 < b < 1.5 (see, e.g.. IVinkll2013 ). The condition 
(|12|l predicts a less efficient plasma heating in the CR-precursor of the faster shocks. 

Heat conduction that we did not discuss here may play a role in the electron 
temperature distribution in the shock precursor (c.f. lBreech et al.ll2009r ). The nonlinear 
dynamics of the th e CR-driven magneti c fiuctuations in the shock precursor deserves 
thorough modeling. iMalkov et al.l ( 20121 ) have obtained fully nonlinear exact solutions 



of the ideaJ ID-MHD supported by the CR return current. The solutions occur as 
localized spikes of circularly polarized Alfven envelopes (solitons or breathers). The 
sufficiently strong solitons in the model run ahead of the main shock and stand in 
the precursor, being supported by the return current. The CR-shock precursor in the 
model is dissipationless. 

The electron and ion temperatures in the shock precursor determine the injection 
of particles into the CR acceleration regime. The temperatures can be tested by optical 
spectroscopy of supernova shocks. CR-modification of shocks with modest speed of a 
few hundred km s~ may yield lower post-shock temperatures and thus make the 
post-shock flow switch to a radiative regime. 



5 Spectroscopy of a CR-modified radiative shock p 

Consider a one-dimensional flow around a coUisionless shock consisting of three zones: 
a) the pre-shock, where the unperturbed interstellar matter is preionized and preheated 
by the radiation (and energetic particles) generated in the downstream and where 
strong fluctuations of magnetic field may be generated by the CR anisotropy, b) the 
thin shock front (a "viscous jump" ) where a substantial part of the kinetic energy of the 
bulk upstream flow is converted into energy of thermal motions, and c) the post-shock. 



^ This section uses some results of numerical code SHELLS currently developed by 
A.M. Bykov, A.E. Vladimirov, and A.M. Krassilchtchikov. This code is used to model the 
steady-state structure and broadband continuum and line emission spectra of radiative shocks. 
The model accounts for the impact of efficient CR acceleration on the shock compression 
ratio and the postshock flow. The code incorporates modern atomic data and allows us to 
consider in detail the non-equilibrium microscopic, thermodynamical and radiative processes 
that determine the plasma flow. 
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where the hot flow cools down, radiating continuum and line emission. We discuss here 
a class of shock flows in partly ionized media, the so-called radiative shocks, where the 
power radiated away from the post-shock flow is a siz eable fraction of the total kinetic 
and magnetic power dissipated at the shock (see, e.g.. lSpitzeiill978l : lDraine and McKeg 
I1993I I. 

Without a signiflcant impact on accuracy, one may assume that the pre-shock is 
isothermal, and the flow is in a steady state (see also a comment about the equilibrium 
between the pre-shock ionization state and the ionizing flux at the end of Section [2]) . 
However, the ionization state of the plasma is non-uniform: the ionization level increases 
toward the subshock, as the gas absorbs the ionizing radiation flux emerging from the 
hot downstream region. This photoionization does not have a signiflcant impact on 
the thermal state of the gas, because the hot electrons produced by photoionization 
do not have sufficient time to coUisionally equilibrate with the atoms and ions. CR 
particles m ay contribute to gas ionization and hea ting in the interstellar clouds and 
shells ('e.g.. lSpitzeilll978l : [Bozhokin and EvkovHigQ-d l. The effect of CRs may be impor- 
ta nt in the vicinity of fas t Balmer-type interstellar shocks as it has been demonstrated 
bv lMorlino et al.l (|2012al ). who did not account for the photo processes, though. Ra- 
diative shocks are expected to occur in relatively dense environment and, therefore, 
have velocities typically well below 1,000 km s~ . In such an environment, photo pro- 
cesses can provide a high ionization degree of the upstream gas for sufficiently high 
shock speeds. If CR acceleration is weak, the photoionization becomes significant for 
shock velocities above 100 km s~ . However, if CR acceleration is efficient, the down- 
stream temperature is reduced and the shock compression is increased, which means 
that strong photoionization in the precursor occurs at much greater shock speeds See 
Fig. U for an illustration of the effect of CR acceleration on precursor ionization. We 
demonstrate below the effect of the CR fiuid on the shock compression and spectra of 
radiative shocks. 

The jump conditions at the shock are either given by the Rankine-Hugoniot equa- 
tions, or, if the shock structure is as sumed to be modified by accelerated particles (e.g., 
lBvkovll2004l : IVladimirov et al.ll2008l . and references therein) , the compression ratio can 
be parameterized by the fraction Qesc of bulk fiow energy carried away by these parti- 
cles. The subshock is the standard gas viscous shock of a Mach number Algub- For that 
simplified two-Huid model of a strong CR-modified shock the effective ion temperature 
in the downstream T. can be estimated for a shock of a given velocity, if rtot and 



''sub are known: 






Tg'-totl^sh) 


where (t){Msuh) = 


27gMsub - (7g - 1) 
(7g-l)-Ms'ub + 2 



(13) 

Single fluid strong shock heating represents the limit A4snh = Als 3> 1, since there 
is no precursor in that case the temperature behind a strong shock is determined by 
the standard scaling 

y(2) ^ 2 . p^^ M«sh = 1-38 • 10^ ^.^8 (K), (14) 

(7g + 1)^ 

In single-fluid systems the compression ratio rtot = ^sub ~^ (7g + l)/(7g ~ 1) does 
not depend on the shock velocity and Ea. (|13|) reduces to Ea. (|14p . However, in multi- 
fluid shocks the total compression ratio depends on the shock velocity and could be 
substantially higher than that in the single-fluid case. Consequently, the post-shock 
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temperature in a multi-fluid shock is lower than the post-shock temperature in a 
single-fluid shock of the same speed. This allows us to determine the CR accelera- 
tion efficiency using observations of post-shock temperatures and shock speeds, or the 



tion emciency using observations or post-snocK temperatures ana snocK speeas, or tne 
entropy proflles in the accretion shocks of clus ters of galaxies ( Bvkovll2005l : iBvkov et alj 



I2OO8I : iBriiggen et al.ll2012l : iFuiita et al.ll2013l ). It is convenient to introduce the scal- 
ing rtot(tJsh) oc V u to describe the different cases of strong shock heating ( Bvkov et alJ 



I2OO8I '). Then from Eq. ^ T.^'^^ ex ?i(Xsub) '"fh^ ^' ■ The subshock Mach number X^ub 
depends, in general, on Als and A4a,. Thus, the index a approximates the velocity de- 
pendence of ^(TVfsub) o^ v^h- Finally, HT. ' oc w^j^, then the index a = 2(1 — ^)-|-cr . For 
the case of sh ock precursor heating by CR generated Alfven waves, the index a ~ 1.25 
(|Bvkovll200a ). The effects of neutrals due to charge exchange in the shock downstream 
with heated ion s that results in a flux of high-veloc ity neutrals that return upstream 
were studied bv lBlasi et al.l (2013) and lOhiral ( 2012h . They found that the return flux 



of neutrals may result in the reduced shock compression ratio and spectral steepening 
of test particles accelerated at the shocks slower than about 3000 km s~ . The return 
flux of neutral atoms may also affect the radiation spectrum of the post-shock flow 
that we are modeling. 

The ratio of the electron to ion temperature immediately after the subshock, 
Te/Ti = 5e, is considered as a free parameter in our model. It can be estimated from ob- 
servations as discussed in Section [S] It is varied in the range from y^nie/mp ~ 0.023 to 
1. Observations typically show low value s of 5e in shocks faster than 1000 km s~ , and 
values clo ser to 1 in slow er shocks (e.g.. iGhavamian et al.ll200ll . 120071 : iRakowski et al.l 
[2OO8; Hol der et al.ll2010l ). 

The post-shock plasma is treated as a stationary two-temperature single fluid flow 
consisting of ions, electrons, and neutral atoms t hat are ideal nonrelativistic gases. The 



neutr al and ion temperatures may actually differ ( Heng et al.ll2007l : lvan Adelsberg et al.l 
[2003), but that is only important just behind a very fast shock in partially neutral gas, 
and it matters mainly for diagnostics based on the Ha line proflle. For the considered 
here ranges of shock speeds and matter densities the temperature equilibration scales 
are below 10 cm~ , while the line emission zone typically lies well above 10 cm~ 
downstream. 

Let z and ^ be the average charge and mass of an ion, (^ = me/i-i - the average 
electron to ion mass ratio. Then 

Ue = zrij, pe — merie = zmeUi, (15) 

Pi = MJij, 

p = ni{zme + m) = nip{zC + 1) = nm/M{z), 

M{z) = l/{zC+l). 

Let po,vo, fi'iid Bq be the values of density, velocity, and frozen-in transverse magnetic 
field just before the shock, and Va,T°' be the values of flow velocity and ion temperature 
in the immediate post-shock defined by jump conditions. 

Then, the flow evolution downstream can be described by the following system of 
equations. 



pv = const — povo, (16) 

pv'^ +Pe+ Pi + Pm + PcR = const = npaVQ, (17) 
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3 dT, dv SrUeUe 

3 dTe dv SmeUe , rj, •, . 3 dz 

r^"^ = -"^^^ dJ + —-,^^^ -Te)-A- ^n^TeV-, (19) 

where pe + Pi = nikB{zTe + T^) — pkgM^zTe + T^)/^, the magnetic pressure pm = 

Blp^/ifi-Kpl), n = «a/«o+Bo/(8'rpo"0«a)+fcsMTf (z4 + l)/(«0"aM)+Qosc/2, PCR = 

(Qcsc/2) ■ poVQ{va/v) ' , yl is the coohng function, Te^ is the electron-ion equiUbration 
time, and the last term of equation H19|) denotes electron cooling due to ionization. 
The cooling term A is calculated as 

A = n'AcoU + rie y^^in^^jat^jElj" + rnjCijhuj) - (20) 

1,0 •'"^ 

where i denotes the chemical element and j denotes the ionization state of an ion (j=0 
corresponds to a neutral atom). Here A^^u is due to electron-ion collisions includ- 
ing electron bremsstrahlung and line emission of the ions excited by electron impact, 
though it does not include emission due to radiative recombination; a,; j is the recom- 
bination rate {j -|- 1} — >■ {j}; E^''^ is the average energy on the recombining electrons; 
Ci , is the rate of coUisional ionization; af - is the photoionization crossection of the 
ion state j; Jy is the angle-averaged density of ionizing radiation at frequency v. 

Ju{x) = -!- / 27r Iy{iJ.,x) dfj,. (21) 

where p = cos{0), 6 is the angle between the normal to the shock front and the direction 
of emitted photons. 

The evolution of the ionization state j of an ion i in the downstream flow is deter- 
mined as 



dn. 



V- 



■ij 



dx 



ne (ni.j^iCij^i - ni^jCi,-j - ni^^a^j + riij^iaij^i) + (22) 

, \ ^ T rion f-irion , -irrec \ , ir^GC 

+ 2^ ris nij_iV/5 J _,■_! -nij(Ks i J -I- Ks i j) +nij+il/s_ij + i 

s=H,He,He+ 

where Ri^j is the photoionization rate, Vs°^ and 14*°" are the rates of ionization and 
recombination via charge exchange reactions with the ion s, Ci_j — Ci_j + Cf^^°, 
where Cf^ ° is the autoionization rate, Oij = a^.j -t-ftj*:, where aj*:- is the dielectronic 
recombination rate. 

To obtain the ionizing radiation field a transfer equation can be solved both in the 
downstream and in the upstream: 

cosS • — -^ — ju — Kijlu, (23) 

dx 

where 9 is the angle between the normal to the shock front and the direction of emitted 
photons. In this equation, the absorption coefficient i^iy ~^- ,ni,j(j^- is determined 
by bound-free transitions in all ion species. The ionizing emission is generated as a) 
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Fig. 4 Hydrogen and helium pre-shock (photo)ionization in the shock upstream for four sets 
of shock parameters simulated with the model of radiative shock accounting for the effect of 
CR escape on the shock compression r (the CR-modified case is shown in the top left panel). 
Within the presented model r = 6.2 corresponds to 30% of the flow energy being converted 
into CRs. 



permitted ultraviolet and optical lines excited by an electronic impact, b) recombina- 
tion line cascades of hydrogen and helium, c) free- free continuum of electrons scattering 
at ions, d) 2-photon continuum emission of H- and He-like ions where metastable levels 
are coUisionaly populated, e) recombination continuum emission. 

The relative contribution of each of the mechanisms to the ionizing photon field 
is illustrated in Figure O The dominating 2-photon continuum emission comes mainly 
from He I and He H and the dominating ultraviolet line emission comes mainly from 
He I, He H, and oxygen ions upto O V. It should be noted that in the presented model 
the electron temperatures just after the shock are assumed to be low (Se ~ 10~ ) and 
equilibrate with the ion temperatures downstream via Coulomb collisions. 

All the ionizing lines except Lya were considered optically thin. The optically 
thick case of Lya was treated a ccording to the standard 2-level formalism adopted 
from lMihalas and Mihalai (| 19841 '). 
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Fig. 5 Production rate of ionizing photons (hu > 1 Ry) in the downstream of a 130 km/s 
shock penetrating into a 10 cm~'^ medium of solar abundance. Here "fb" denotes recombination 
continuum emission, "fF' denotes free- free continuum of electrons scattering at ions, "2ph" 
stands for 2-photon continuum emission of H- and He-like ions, and "line" denotes permitted 
ultraviolet lines excited by an electronic impact. The dashed line illustrates the simulated 
electron temperature profile in the downstream. The photons produced in the optically thin 
part of the downstream are ionizing the upstream flow. 



Figures U through [8] and Tables [T] and [2] demonstrate the results of modeling with 
SHELLS. The impact of efficient CR acceleration on the shock is accounted for in all 
calculations by using a reduced downstream temperature and the CR pressure in the 
shock downstream. 

The steady-state pre-shock ionization structure is shown in Figure [4] for shocks 
with speeds of 100, 130, 200, and 310 km s^ illustrating also an effect of the CR fluid 
on the pre-shock ionization. The plasma temperature and gas compression just behind 
the CR modified shock affect the ionization structure. One can see in Figure[4]that the 
ionization structure of 310 km s~ shock modified by CRs may be similar to that of 
100 km s^ shock without CRs fiuid effect. This Figure shows the ionization level of 
hydrogen and helium; however, all abundant chemical species up to Fe are tracked in 
the calculation. At speeds of approximately 130 km/s and above, hydrogen upstream 
becomes completely ionized, and the extent of the radiative precursor is determined 
by the competition between hydrogen and helium recombination and photoionization. 

Figure E] illustrates the downstream cooling region of a 200 km/s shock. Three lines 
correspond to three models with different assumptions regarding the efficiency of CR 
acceleration. For efficient acceleration, energy conservation requires greater compres- 
sion ratio and lower downstream temperature. The three cases shown here illustrate 
shocks with compression ratios of 4.6, 5.3 and 6.2. 
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Fig. 6 Gas temperature (solid lines) and density (dotted lines) in the downstream cooling 
region of three 200 km/s shocks. Three lines correspond to three models with different as- 
sumptions regarding the efficiency of CR acceleration. Within the presented model r = 6.2 
corresponds to Qosc = 30% of the fiow energy being converted into CRs, r = 5.3 corresponds 
to Qgsc = 20%, and r = 4.6 corresponds to Qesc = 10%. 



Shown in Figure [7] is the result of a calculation with shock and ISM parameters 
plausible for the region of interaction of the north-eastern shell of the SNR IC 443 with 
an interstellar cloud. Unperturbed gas density is uq = 15 cm" , shock velocity vq = 
130 km/s, and perpendicular magnetic field Bg = 3 fj,G are assumed. The temperature 
and total density of the plasma are plotted as a function of gas column density from 
the subshock. These quantities are then used to calculate the radiative transfer of 
emission lines in the cooling flow, which can be used for shock diagnostics. Figure |5] 
shows the profiles of two of the most prominent infrared lines (with respect to estimated 
backgrounds) for this case: C II 157.7 fim and N II 205.3 fim. 

The model presented in this Section may be used to diagnose various parameters 
of radiative shocks with observed infrared lines, including the cosmic ray acceleration 
efficiency. 
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Fig. 7 Post-shock flow profile for shock and ISM parameters plausible for the region of inter- 
action of the north-eastern shell of IC 443 SNR with an interstellar cloud. Here the upstream 
density no = 15 cm"'', the shock velocity vo = 130 km/s, and the perpendicular magnetic 
field Bo = 3 fiG. 



6 Formation of CR spectral breaks in partly ionized shock precursors 



6.1 Gamma-ray observations 



The neutral component of the interstellar medium may manifest itself in the gamma- 
ray emission of some coUisionless shocks. 

The recent Fermi-LAT observations of the so-called molecular SNRs W44 and 
IC 443 (|Abdo et al.ll2010all2010bl : lAckermann et al.ll2013l ') indicate that the spectra of 
the gamma-ray producing protons (integrated over the emission region) are typically 
steeper than the DSA predictions for the spectra of the CRs confined in the acceleration 
region. The steep photon spectra has been found in the high energy gamma-ray spectra 
of some other remnants measu red by, e.g., t he CANGAROO (Enomoto ct al. 20021. 
H.E.S.S JAharonian et aLlbOOd ) and MAGIC (|Carmonall2oTil ) atmospheric Cherenkov 
telescopes. 
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Fig. 8 [C II] 157.7 /^m and [N II] 205.3 irca lino profiles for the backgrounds of 140 MJy/sr and 
130 MJy/sr, respectively, estimated for the north-eastern shell of IC 443 SNR with Herschel 
Spot (for the date 22/09/2010). 



Particles accelerated by DSA mechanism generally comprised two distinct popu- 
lations - CRs confined in the accelerator and the CRs escaping the system and these 
two populations have very different spectral s hapes. The observed ga mma-ray emis- 
sion is a sum of the two contributions (see, e.g.. lEllison and Bvkovll2011h . The common 
feature of the molecular SNRs is a significant amo unt of dense molecu lar gas in their 
surroundings. It has been argued on this ground ( Malkov et al.ll2005l ). that when a 
SNR interacts with a dense molecular cloud complex, the conditions for particle con- 
finement to the shock are different from those adopted in conventional DSA modeling. 
Since the propagation of resonant Alfven waves is inhibited by ion-neutr al collisions, 
some particles are not confined and so escap e the emission volume (e.g. iDrurv et al.l 
Il996l : iBvkov et aLllJOOol : iMalkov et al.ll2005l '). The accelerated particle partial escape 
should result in a spectral break in their spectrum and thus, in that of the observed 
gamma emission. In a clear-cut case of such a limited CR confinement, the spectral 
index at the break should change by exactly one power Aq = 1 due to an effective 
reduction of the particle momentum space dimensionality by one, since particles are 
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confined in coordinate space only wlien tliey are witliin a slab in momentum space 
oriented perpendicular to the local magnetic field. 

The most convincing evidenc e for the breaks of index one are the recent Fermi-LAT 
and A GILE observati ons of W44 (lAbdo et al.ll2010al : lu"chivama et al.ll2010l : [Giuliani et all 



12011 



) (re-analyzed in dMalkov et al.lboill 'l'l and MAGIC observations of the SNR W51C 



m 



( Carmonall201ll : lAleksic et al.ll2012l ) These observations are encouraging in that they 



unambiguously confirm the breaks. They demonstrate departure from the traditional 
DSA models in fully ionized plasma where the proton spectrum is a single power law 
with an exponential cutoff. A possible explanation to the spectral break in the gamma- 
ray emission of these objects is that neutral atoms modify the nature of wave-particle 
interactions, leading to a spectral break. This idea is discussed below. 



6.2 The DSA mechanism 

A great deal of the success of the diffusive shock acceleration mechanism (DSA) is 
due to its ability to reproduce the observed power law spectra associated with possi- 
ble CR accelerators, such as SNR shocks. For a classical step-like shock, propagating 
into an ionized medium that supports sufficient CR scattering on both sides of the 
shock, the mechanism predicts an almost universal power-law energy spectrum for the 
accelerated CRs, ex E~'^ , with an index q = (r -I- 2) / (r — 1). It thus depends only 
on the shock compression r, that is, however, close to four for strong shocks. The 
backreaction of the accelerated CRs on the shock structure may change this result 
noticeably, but not dramatically. For strong shocks, 7V( S> 1, the spectrum hardens 
from g = 2tog~1.5, at most. This change is largely due to an enhanced com- 
pression of the shock that results from the CR escape flux and the reduction of the 
adiabatic index of the CR/thermal-plasma mixture compared to the ordinary plasma 
(i.e., 7 — 5/3 — >■ 4/3). More problematic is to soften the CR spectrum by their back- 
reaction effects on the shock environment. This is intuitively understandable, since 
spectrum steepening should diminish the backreaction coming from the softening of 

the equation of state (7 — >■ 4/3) and, especially, from the CR escape flux. 

The recent Fermi-LAT observations of the SNRs W44 and IC 443 (|Abdo et all 



I2010all2010bl l also urge a second look at the DSA mechanism (see Section [6.1 [I 



6.3 Mechanism for a spectral break 

When a SNR shock approaches a molecular cloud (MC) or a pre-supernova swept-up 
shell with a signiflcant amount of neutrals, confinement of accelerated particles deteri- 
orates. Due to the particle interaction with magnet ic fluctuations, confinement requires 
their scales to be similar to the particle gyroradii ( Drurvlll983l : lBlandford and Eichleii 



119871 ). However, strong ion-neutral collisions make the wave-particle interactions more 
sensitive to the particle pitch angle, which can be understood from the following con- 
sideration. 

While the waves are in a strongly ionized (e.g., closer to the shock) medium they 
propagate freely in a broad frequency range at the Alfven speed Va = -B/V^TrPi with 
the frequencies to = kV^- Here k is the wave number (assumed parallel to the local 
field B) and Pi is the ion mass density. As long as the Alfven wave frequency is higher 
than the ion- neutral collision frequency j^j„, the waves are only weakly damped. When, 
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on the other hand, the ion-neutral colhsion frequency is higher (deeper into the cloud) , 
neutrals are entrained by the oscillating plasma and the Alfven waves are also able to 
propagate, albeit with a factor ypi/fKi < 1 lower speed, where po is the neutral density. 
The propagation speed reduction occurs because every ion is now "loaded" with po/Pi 
neutrals. Now, between these two regimes Alfven waves are heavily damped and even 
disappear altogether for sufficiently small Pi/po '^ 0.1. The evanescence wave number 
range is then bounded by fci — i' i„/2VA and fc2 = 2^/|Oj/po t ^i„,/V/^. These phenomen a 



have been studied in detail in ( Kulsrud and Pearcd 196E : IZweibel and Shulll 
and specifically in the context of th e DSA in ( Volk et al.Ml98ll : IPrurv et al 



lisi) 



19961 : 



iBvkov et al.ll200Cl : iReville et al.l 120081 1. Now we turn to their impact on the particle 
confinement and emissivity. 

In the framework of a quasilinear wave-particle interaction, the wave number k is 
approximately related to the parallel (to the magnetic field) component of the particle 
momentum pii by the cyclotron resonance condition kpu/m = ±iJc, where the (non- 
relativistic) gyro-frequency ujc = eB/mc. Note that the appearance of pji = p/x, where 
p is the cosine of the pitch angle instead of the often used "sharpened" (|Skillindll975l 'l 
resonance condition kp/m — ±ajc is absolutely critical for the break mechanism. The 
frequency range where the waves cannot propagate may be conveniently translated into 
the parallel momentum range 

Pi < |P||| <P2, (24) 

with 



Pi = 2VAmu}c/ui„, P2 = -j\/po/Pi >Pl- (25) 

That a spectral break must form at the photon energy corresponding to the particle 
momentum p = pi = Phr, can be readily explained as follows. The 'dead zones' pi < 
PlI < P2 imply that particles with pii > pi do not turn around (while moving along 
the magnetic field) and escape from the region of CR-dense gas collisions at a pu/p 
fraction of the speed of light. More specifically, particles with pi < Ipu < p2 escape 
because they are not scattered, whereas particles with pii > P2, because they maintain 
the sign of pu , even though they scatter but cannot jump over the gap pi < pii < P2. 
An exception to this rule are particles with sufficiently large pj_ that can be mirrored 
across the gap or overcome it via the resonance broadening. 

The break can also be explained in terms of the confinement times of different 
groups of particles introduced above if we assume a low density pre-shock medium 
with clumps of dense, partially neutral, gas. Particles with pii > pi spend only short 
time Tosc ~ Lc/c (where Lc is the size of the clump) inside the gas clumps. They 
propagate ballistically and their scattering time is assumed to be infinite, as there are 
no waves they can interact with resonantly (pi < pii < P2) or they cannot change 
their propagation direction ( pii > P2)- Particles with pJ < pi are, on the contrary, 
scattered intensively in pitch angle, they frequently change their direction, and so sit 
in the clump for r^ouf ^ L(./n '-^ L^/c Tsc- Here Tsc is their pitch-angle scattering time 
and K is the associated diffusion coefficient. Not only r^-^^f 3> Tcsc is required, i.e., 
Tsc ^ Lc/c, but also Tconf > Lc/Ush, which means that the shock precursor is shorter 
than the clump n/Ush ^ icR < ^c (here C/gh is the shock velocity, and icR is the 
thickness of the CR front near the shock). The last condition ensures that particles 
with p|| > pi that escape through the clump after having entered it from the shock 
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Clumps of dense gas 




SNR-shock 



Fig. 9 SNR shock propagating into dense gas environment. The filling factor of the 
gas clumps is small, while some of thorn may be larger than the thickness of the CR layer near 
the shock front. 



side, will not interact with the shock after they exit through the opposite side of the 
clump, thus escaping upstream, Fig[§l The reason for that is a low level of Alfven wave 
turbulence ahead of the CR precursor. We also assume that the ambient magnetic field 
does not deviate strongly from the shock normal, in order to allow these particles to 
escape through the far side of the clump. 

While particles with p > pi escape from the regions of enhanced gamma radiation 
(high gas density), an initially isotropic distribution of accelerated particles is main- 
tained only in a slab in momentum space \pu < pi and becomes thus highly anisotropic 
(a 'pancake' distribution). What matters for the integral emission, however, is a locally 
isotropic component / of this new proton distribution. It can be introduced by re- 
averaging the 'pancake' {\pu < pi) distribution in pitch angle, / (p) = J f (p, fi) dfi, 
and is readily obtained assuming that particles remaining in the dense gas (those with 
\p\\ < Pi) maintain the flat pitch-angle distribution, i.e. 



Ml 



fip) 



/o (p) d^J' - 



(pi/p)/o (p), P> Pi 
h (p) , p < Pi 



(26) 
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where /o (p) is the initial (isotropic) distribution function and m — min{pi/p, 1}. 
Thus, the slope of the particle momentum distribution becomes steeper by exactly one 
power above p — pi = p]^,.. In particular, any power-law distribution ex p~'^ , upon 
entering an MC, turns into p~'~ at p > pjj^, and preserves its form at p < p^^.. 

Note that the broken power-law spectrum can only be maintained if the filling factor 
/gas of the dense gas with the significant wave evanescence interval (pi , P2) is relatively 
small, /gas ^ 1, so that the overall particle confinement and thus the acceleration are 
not strongly affected. If, on the contrary, /gas ~ 1, the resonant particles would leak into 
the (pi,P2) gap and escape from the accelerator in large amounts, thus suppressing the 
acceleration. Note that particles with sufficiently high momenta p > P2Bq/5B, where 
SB/Bq is the effective mirror ratio of magnetic perturbations, can "jump" over the 
gap. The primary p~'^ slope shoul d then be restored for such particles . Recent MAGIC 
observations of the SNR W51C jCarmonal [2OIII : lAIeksic et al.lboid ') indeed point at 



such spectrum recovery at higher energies. It should also be noted, that the Aq = 1 
break index is a limiting case of identical gas clumps. The integrated emission from an 
ensemble of clumps with different pi and p2 should result in a more complex spectrum. 



6.4 Break momentum 

While the one power spectral break in the pitch-angle averaged particle distribution 
seems to be a robust environmental signature of a weakly ionized medium into which the 
accelerated particles propagate, the break momentum remains uncertain. According to 
eQ. H25|) . Pf,,- (= Pi) depends on the magnetic field strength and ion density as well as on 
the ion-neutral collision rate Uin ~ no (crV). Here (aV) is the product of the collision 
cross-section and collision velocity averaged over the thermal distribution. Using an 
approximation of ( Draine and McKed 119931 : iDrurv et al.lll99(i ) for (cV), pi„, can be 



estimated as 



PbJmc ~ WB^T^ °-\ ^n^ ^''^ . (27) 



Here the gas temperature T4 is measured in the units of 10 A", magnetic field B^ -in 
microgauss, np a-nd n^ (number densities corresponding to the neutral/ion mass den- 
sities po s^nd Pi) -in cm~ . Note that the numerical coefficient in the last expression 
may vary depe nding on the average ion and neutral m asses and can be higher by a 
factor of a few ( Kulsrud and Pearcdll969l : lNakandll984l ) than the estimate in ea. H27|l 



suggests. The remaining quantities in the last formula are also known too poorly to 
make an accurate independent prediction of the position of the break in the gamma ray 
emission region. Those are the region s near the blast wave where complicated physica l 



processes unfold, as discussed earlier ( ShuU and McKeal 19791 : IPraine and McKeal 19931 : 



iBvkov et al.ll200Cl ). Also important may be the ionization by the low energy CRs accel- 
erated at the blast wave. However, as their diffusion length is shorter than that of the 
particles with p > p;,^, we may assume that they do not reach the MC. Pre-ionization 
by the UV photons can also be ignored for th e column density A'^ > 10 cm~ ahead 
of the shock beyond which they are absorbed |Uchiyama et al.ll20ld ). IUchivama et al.l 
(2010|), using the earlier data from ( Reach et al.ll20Q5l ) have also analyzed the parame- 



ters involved in eq.((27| and found the above estimate of pf,r to be in a good agreement 
with the spectral break position measured by the Fermi-L AT. Nevertheless, w e may 
run the argument in reverse and use the Fermi observations ( Abdo et al.ll2010al ) of the 



gamma-ray spectrum of SNR W44 to determine the break momentum in the parent 
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particle spectrum and constrain tlie parameters in eq.([27|. Since we also know the 
amount of the slope variation Aq, we can calculate the full spectrum up to the cut-off 
energy. 

It should also be noted that in reality the break at p = p]-,r is not infinitely sharp for 
the following reasons. The break momentum may chan ge in space due to variations o f 
the gas parameters (eq.[2Z]), the resonance broadening (|Dupreelll966l : lAchterbereill98ll ) 



near p = pi = Pbr (so that particles with p ~ pi are still scattered, albeit weakly) and 
other factors, such as the contribution of small gas clumps with Lc <^ ^CRi FigEl The 
small clumps are submerged in the CR front and the CRs that escape from them are 
readily replenished. Note that this effect may decrease the break index Aq. However, 
the conversion of the parent proton spectrum into the observable gamma emission 
introduces a significant smoothing of the break, so that even a sharply broken proton 
spectrum produces a smooth gamma spectrum. It provides an excellent fit to the Fermi 
gamma-ray data wi thout an ad hoc proton break smoothing adopted by the Fermi-team 
( Abdo et allboiOal ) to fit the data. 



7 Particle and gamma-photon spectra in molecular cloud SNRs 

To calculate the particle spectra, we need to determine the degree of nonlinear mod- 
ification of the shock structure. In principle, it can be calculated consistently, given 
the shock parameters and the particle maximum momentum, pmax- In the case of a 
broken spectrum, pf,^ likely plays the role oi pmax, as a momentum where the domi- 
nant contribution to the pressure of accelerated particles comes from, thus setting the 
scale of the modified shock precursor. Note that in the conventional nonlinear (NL) 
acceleration theory, the cut-off momentum pmax plays this role, because the nonlinear 
spectra are sufficiently flat so as to make the pressure diverge with momentum, unlike 
the broken spectra. 

One of the best documented gamma emission spectra comes from the SNR W44, 
so we flt these data using the above mechanism of the spectral break. The break 
in the photon spectrum is observed at abou t 2 GeV, which pl aces the break in the 
proton distribution at about pi,,. ~ IGeV/c I Abdo et al.ll2010al ). For the strength of 



the break Aq = 1, the spectrum above it is clearly pressure converging, so that the 
shock structure and the spectrum may be calculated using this break momentum as 
the point of the maximum in the CR partial pressure. More specifically, once the break 
momentum is set, one can use an analytic approach (.Malkov and Drurv 200l|) for a 
stationary nonlinear acceleration problem using pt,^ as an input parameter. 

Apart from p^r-, the nonlinear solution depends on a number of other parameters, 
such as the injection rate of thermal particles into acceleration, Mach number, the pre- 
cursor heating rate and the shock velocity Vs . Of these parameters the latter is known 
reasonably well, Vs ~ 300 km/s, the injection rate can be either calculated analytically 
for the parallel shock geometry JMalkov and Drurv^l2001^ . or inferred from the simu- 
lations ( Gargate and Spitkovskvll2012h . while the other parameters are still difficult to 



ascertain. Fortunately, in sufficiently strong shocks the solution either stays close to 
the test particle (TP) solution (leaving the shock structure only weakly modified) or 
else it transitions to a strongly modified NL-solution regime. The TP regime typically 
establishes in cases of moderate Mach numbers, low injection rates and low pmax (now 
probably closer to Pf,,.), while the NL regime is unavoidable in the opposite part of the 
parameter space. 
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In the TP regime the spectrum is close to a power-law with the spectral index 2 
throughout the supra-thermal energy range. In the NL regime, however, the spectrum 
develops a concave form, starting from a softer spectrum at the injection energy, with 
the index q ~ [rs + 2)/{rs — 1) > 2, where r^ < 4 is the sub-shock compression ratio. 
Then it hardens, primarily in the region p ~ mc, where both the partial pressure 
and diffusivity of protons change their momentum dependence. The slope reaches its 
minimum at the cut-off (break) energy, which, depending on the degree of nonlinearity, 
can be as low as 1.5 or even somewhat lower if the cut-off is abrupt. The question now 
is into which of these two categories the W44 spectrum falls? Generally, in cases of 
low maximum (or, equivalently, low spectral break pt.^ < 10) momentum, the shock 
modification is weak, so the spectrum is more likely to be in a slightly nonlinear, almost 
TP regime. On the other hand, there is a putative indication from the electron radio 
emission that their spectrum may be close to qe ~ 1.75, which could be the signature of 
a moderately nonlinear acceleration process. It should be remembered, however, that 
this is a global index across the W44 remnant. There are resolved bright filaments where 
a canonical a = —0.5 spectrum, corresponding precisel y to the TP parent electron 
spectrum with qe — 2, is observed (|Castelletti et al.ll2007h . Moreover, there are regions 
with the positive indices a < 0.4 which cannot be indicative of a DSA process without 
corrections for subsequent spectral transformations such as an absorption by thermal 
electrons. These regions may contribute to the overall spectral hardening discussed 
above, thus mimicking the acceleration nonlinearit y. Finally, se condary electrons give 



rise to the fiattening of the radio spectrum as well (|Uchivama c t al. 201Q) 



The above considerations somewhat weaken the radio data as a probe for the slope 
of the electron and (more importantly) for the proton spectrum. Therefore, the exact 
degree of nonlinearity of the acceleration remains unknown and one can consider both 
the TP and weakly N L regimes in ca l culations of the photon spectra, generated in p — p 
collisions. Specificallv. iMalkov et al.l ( 20051 . 1201 J l calculate the tt production rate and 
th e gamma-ray emissivity. In so doing, they adop t numerical recipe described in detail 
in ( Kamae et al.ll2006l : lKarlsson and KamaduOOa ) . The physical processes behind these 



calculations are (i) collisions of accelerated protons with the protons of the ambient 
gas resulting in pp — >■ it reaction (ii) decay of tt - mesons to generate an observable 
gamma emission spectrum. 

An example of such calculations is shown in Fig llOl The best fit to the Fermi and 
AGILE data is provided by a TP energy distribution (oc E~ ) below p^^, ~ 7GeV/c 
with the spectrum steepening by exactly one power above it. The spectrum steepening 
is perfectly consistent with the proton partial escape described above (with no parame- 
ters involved). For comparison, a weakly NL spectrum (guided by the inferred electron 
spectrum with qe « 1.75, is also used for these calculations (dashed line in Fig llOp . but 
its fit would require a somewhat stronger break {Aq > 1) or a low momentum cut-off, 
i.e. at least one additional free parameter. It is seen that the mechanism for a break in 



the spectrum of shock accelerated proton s suggested in IIMalkov et al. 2005 ) provides a 
good fit to the recent (|Abdo et aLlboiOal ) Fermi-LAT and AGILE (|Giuliani etaHl201ll l 



observations of the SNR W44. Of course, in assessing consistency of the suggested spec- 
tral break phenomenon with the observed spectrum, the errors in the data must be 
taken into account. The vertical error bars near the break at 2 GeV, are fairly small 
(comparable to the size of the symbols used to represent the data in Fig. llOp . More 
significant appe ars to be the energy dispe rsion. However, in the most recent Fermi-LAT 






publication f see I Ackermann et al.l ( 20131 ) including the Supplementary Online Materi- 
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Fig. 10 Gamma radiation spectra. Photon spectra resulting from tt" decay and calcu- 
lated for two different parent proton spectra compared against the Fermi (circles) and AGILE 
(squares) data. Solid line: a test particle acceleration regime with the spectral index q = 2 
below the break and q = 3 above the break at pj^ = 7GeV/c. Dashed line: a moderately 
nonlinear acceleration regime corresponding to the spectrum with g ~ 1.75 and q ~ 2.75 below 
and above the break, respectively. Cut-offs are placed at 3 00 GeV for TP- and 100 GeV, for 
NL-s pectrum. Fermi and AGILE data are adopted from llAbdo et alj|2010al : iGiuliani et al.l 
120111) , respectively. Both curves are well within the error bars of Fermi and AGILE (not shown 
for clarity), which, in turn, overlap iGiuliani et aljr201ll) . 
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als) the energy dispersion is estimated to be less tlian 5% for these energies, so tliat 
the broken power law is indeed consistent with the data. 

Generally, spectral breaks offer a possible resolution to the well known problem 
that some nonlinear DSA models produce spectra which are considerably harder than 
a simple test particle spectrum, and these are not consistent with the gamma-ray ob- 
servations of some of supernova remnants. However, the nonlinear spectrum - i.e., 
diverging in energy- exhausts the shock energy available for the acceleration as the 
cut-off momentum grows, so that a broken spectrum should form. Broken spectra are 
commonly observed and the old paradigm of a single power-law with an exponential 
upjjcr cut-off is maladapted to the recent, greatly improved observations ( Abdo et al.l 
[2010bi .201o3). Note, that the spectrum of the RX J1713. 7-3946 (jAharonian et al..,2006l l 
is also formally consistent with the environmental break mechanism presumably oper- 
ating in W44 surroundin g but with a higher pi,^ ~ 10 GeV/c and thus with stronger 
acceleration nonlinearity ( Malkov et al.ll2005l ). However, this remnant expands into a 



rather complicated environmen t, so it is difficult to make t he case for hadronic ori 
gin of the gamma- ray emission ( Aharonian et al.ll2004l . 120061 : iKatz and Waxmanll2008l : 



Ellison et al.l 120121 '). The important role of the W44 remnant for the problem of CR 
origin is that this particular remnant seems to be unlikely domin ated by the lepton 



emission due to Brem sstrahlung and inverse Compton scattering ( Abdo et al.ll2010al : 
lUchivama et al.ll201Cl l thus favoring the hadronic origin of the gamma emission and 
bolstering the case for the SNR origin of galactic CRs. 



8 Summary 

CoUisionless shocks are ubiquitous in astrophysical objects and are observed at all 
scales starting from the heliosphere and up to cosmological scale shocks observed in 
clusters of galaxies. Theoretical modeling of these shocks is challenging because shock 
relaxation process involves collective plasma oscillations producing long-lived highly 
nonequilibrium components. Straightforward numerical simulations are an extremely 
resource demanding task because of the very wide dynamical range of scales and time 
required to the resolve all of the long-lived components. Such task is even more difficult 
if neutral atoms and molecules are significant in the partially ionized media. 

In this brief review we have discussed some observational appearance of collision- 
less shocks in partially ionized plasma, described the most important physical processes 
operating in the shocks and outlined the observational perspective of nonthermal com- 
ponents diagnostics of astrophysical coUisionless shocks via multiwavelength observa- 
tions. 

Nonthermal components, i.e., energetic charged and neutral particles and fluctuat- 
ing magnetic fields, can drastically modify the structure of the shock upstream provid- 
ing both deceleration of the plasma flow and also efficient heating of ions and electrons. 
We argue that the processes of turbulence amplification and damping observed in the 
heliosphere may help to understand the microphysics of ion and electron heating by 
cosmic ray driven turbulence in the upstream regions of large scale coUisionless shocks 
observed in the Galaxy and in clusters of galaxies. In turn, radiative signatures of the 
astrophysical shocks can be a unique way to study microscopic phenomena that can 
not be studied in the laboratory plasma on the Earth. 

Ha line diagnostics of Balmer shocks provides valuable information on charge- 
exchange processes of the neutrals with both thermal and non-maxwellian plasma 
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components in the shock downstream. Nonthermal components may reduce the plasma 
temperature comparing to what is expected from the standard single-fluid Rankine- 
Hugoniot jump conditions for a particular shock velocity. In the case of radiative shocks 
this reduction would substantially modify the emission line spectrum coming from 
shock downstream and thus may serve as a valuable diagnostic tool for fast shocks 
interacting with clouds. The effect of neutrals on the MHD wave damping in the 
upstream of supernova shocks interacting with a molecular cloud may explain some 
spectral features in GeV energy regime, recently revealed by Fermi observations. 
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